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Abstract. Given a right factor and a left factor of a Linear Partial Differential 
Operator (LPDO), under which conditions we can refine these two- factor 
factorizations into one three-factor factorization? This problem is solved for 
LPDOs of arbitrary order and number of variables. A more general result for 
the incomplete factorizations of LPDOs is proved as well. 
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1. Introduction 

The factorization of Linear Partial Differential Operators (LPDOs) is an essen- 
tial part of recent algorithms for the exact solution for Linear Partial Differential 
Equations (LPDEs). Examples of such algorithms include numerous generaliza- 
tions and modifications of the 18th-century Laplace Transformations Method [531 
[21 ini HI HH H 120], the Loewy decomposition method [m [H [T3] , and others. 

The problem of constructing a general factorization algorithm for an LPDO 
is still an open problem, although several important contributions have been made 
in recent decades, and different approaches have been applied (see [TTl [H [HI [121 
[2S1 [ini [3 ISl [7] and many others) . Many of the recent approaches are concerned, 
in particular, with explaining the non-uniqueness of factorization: (irreducible) 
factors and the number of factors are not necessarily the same for two different 
factorizations of the same operator. This is commonly illustrated by the famous 
example of Landau |S] , 
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Example (Landau). 

L = (d, + 1+ \ ] o ( + 1 ^— ) o {D, + xDy) = 

= {D^^+xD^y + D^ + {2 + x)Dy)o{D,^ + l) , 

where the second-order factor in the second factorization is hyperboUc and is 
irreducible. 

On the other hand, for some classes of LPDOs factorization is unique. For 
example, there is no more than one factorization that extends a factorization of 
the principal symbol of the operator into co-prime factors [11 . Algebraic theories 
have been introduced to explain this phenomenon theoretically; see Tsarev [21j . 
Grigoriev and Schwarz 13 and most recently Cassidy and Singer [16] . 

Some important methods of exact integration, for example, the Loewy de- 
composition methods mentioned above, require LPDOs to have a number of differ- 
ent factorizations of certain types. Also completely reducible LPDOs introduced 
in TT, which become significant as the solution space of a completely reducible 
LPDO coincides with the sum of those of its irreducible right factors may require 
a number of right factors. 

In earlier work [17] we have exhaustively studied families of factorizations 
for operators up to order 4, and described when the same operator has multiple 
factorizations of the same factorization type, to be more specific, when there exist 
an infinite number of factorizations of the same factorization type, meaning having 
the same symbols of the factors. The first non-trivial example of such families of 
order 4 has been found: 

Example. [TJ The following is a fourth-order irreducible family of factorizations: 

D^^yy = (d^+ " n )iDy+ , , r) iDxy —^—-r{D^ + aDy)], 

V y + ax + pJ\ y + ax + /3J\ y + ax + f3 "I 

where q;,/3 ^ 0. Note that the first two factors commute. So the operator D^xyy 
has a family of factorizations, and every factorization of the family is of the same 
factorization type {X){Y){XY), that is the highest order terms in the first, the 
second and the third factors are D^, Dy and D^y correspondingly. 

In recent work |10| non-uniqueness of a different kind is addressed. There, we 
considered factorizations of different factorization types, and by using invariants 
proved that a third-order bivariate operator L has a first-order left factor of the 
symbol 5*1 and a first-order right factor of the symbol 5*2, where gcd(5i,52) = 1 
if and only if it has a complete factorization of the type (S'i)(T)(S'2), where T — 
Sym(L)/(5i52). Further investigations in the same paper show that a third-order 
bivariate operator L has a first-order left factor Fi and a first-order right factor 
F2 with gcd(Sym(Fi), Sym(i^2)) = 1 if and only if L has a factorization into three 
factors, the left one of which is exactly Fi and the right one is exactly i^2- 
Example. 10 The existence of two factorizations for an LPDO, 

(D^ +x)o {D^y + yD^ + y^Dy +y'^)^A^ {D^^ + {x + y^)D^ + xy^) o {Dy + y) 
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implies the existence of the "complete" factorization of A, 
A={D^ + x)o {D^ + 2/2) o {Dy + y) . 

On the other hand, if the condition gcd(Sym(i^i), Sym(F2)) — 1 fails, then it 
can happen that the "complete" factorization does not exist. 
Example. [lOj 

{D^Dy + 1) o (D^ + 1) = {D,, + 1) o {D^Dy + 1) , 
while Dj.Dy + 1 has no factorization at all. 

In the present paper we have generalized the result of [lO' to the case of 
LPDOs of arbitrary order and of arbitrary dimension. Moreover, a more general 
statement has been formulated and proved for incomplete factorizations of LPDOs. 
We describe the results in terms of common obstacles, which we have introduced 
in [H]. 



2. Preliminaries 

Consider a field K of characteristic zero with commuting derivations di, . . . , dn, 
and the corresponding non-commutative ring of linear partial differential operators 
(LPDOs) K[D] = K[Di, . . . where Di corresponds to the derivation di for 

alH G {1, . . . , n}. In K[D] the variables Z)„ commute with each other, but 

not with elements of K. We write multiplication in K[D] as "o"; i.e. Li o L2 for 
Li,L2 G K[D]. Any operator L ^0 € K[D] has the form 

d 

L=Y, ajD' ,ajeK , (2.1) 
|j|=o 

where J — (ji,...,j„) is a multi-index in N", |J| = ji + •■• + jn, and where 
D"^ = Dl^ . . -Di^ . Further, there exists some J, with |J| = c?, such that aj ^ 0. 
Then d is the order of L. For the case L = we define the order as —00. 

When considering the bivariate case n = 2, we use the following formal nota- 
tions: di = dx, 82 = dy, di{f) = fx, d2{f) = fy, where f G K, and correspondingly 
Di = Dx, and D2 = Dy for ease of notation. 

For an operator L 7^ of the form (|2.ip the homogeneous commutative 
polynomial 

Sym(L) = ^ ajX' (2.2) 

in formal variables X\, . . . ,Xn is called the (principal) symbol, and if X = 0, 
the symbol is defined to be zero. Vice versa, given a homogeneous commutative 
polynomial S G K[X] in the form (12. 2p . we define the operator S G K[D] as the 
result of substituting Di for each variable Xi. 
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3. Main Result 

Since for two LPDOs Li, L2 G K[D] we have Sym(Li o L2) — Sym(Li) • Sym(L2), 
any factorization of an LPDO extends some factorization of its symbol. In general, 
if L G K[D] and Sym(i) ^ Si . . . Sk, let us say that the factorization 

L^Fio-.-oFk , Sym(F,) = S,, Vz £ {1, . . . , fc} , 

is of the factorization type (Si) . . . (Sk)- 

For the second-order hyperbolic LPDOs, which have normalized form 

L = D^Dy + aD^ + bDy + c , (3.1) 

where a,b,c d K, it is common to consider their incomplete factorizations: 

L ^ {D^ + b)o {Dy + a) + h= {Dy + a) o {D^ + b) + k , 

where h =^ c — — ah and k = c — by — ab are invariants of p.ip with respect 
to gauge transformations, L —> g^^Lg, g 0, g G K and are called the Laplace 
invariants. This is an element in the foundation of the classical Laplace-Darboux- 
Transformations Method [3]. 

In [IHllin] a generalization of this idea is suggested. Thus, for A G K[D] with 
Sym(^) = Si . ..Sk, we call [TSl [TO] an LPDO R € K[D] a common obstacle to 
factorization of the type (iS'i)(S'2) . . . (5*^) if there exists a factorization of this type 
for the operator A — R, and R has minimal possible order. 

The following example demonstrates different possibilities for common ob- 
stacles and incomplete factorizations. 

Example. Consider the LPDO 

A4 - DlDl + D^ + Dy + 1. (3.2) 

1. Unique common obstacle and unique incomplete factorization. Consider 
factorizations of A4 of the factorization type {X'^){Y'^). Assume that the order 
of common obstacles is one or less (if we come to a contradiction, we have to 
search then for higher-order common obstacles), and search for common obsta- 
cles in the form Ri — piD^ + qiDy + ri, where pi,gi,ri G K. Thus, for some 
^10, ^01, ^00, /lo, /oi, /oo G K we have 

A4 - {Dl + hoD, + loiDy + loo) o [Dl + fioD, + foiDy + foo) + Ri ■ 

Comparing the corresponding coefficients we have ^lo = ^oi = ^00 = /lo = /oi = 
foo = 0, pi = qi = ri = 1, that is, there is a unique common obstacle and a unique 
incomplete factorization of the factorization type {X'^)(Y'^), 

Ai = DloDl + D^ + Dy + l . 

2. Infinitely many common obstacles and incomplete factorizations. Consider 
factorizations of A4 of the factorization type {X){XY'^). Again assume that the 
order of common obstacles is one or less, and search for common obstacles in the 
form i?2 = PiDx + q-iDy + r2, where P2, 92, 7'2 G K. Thus, for some tooo, gij G K 
we have A4 = {D^ + moo) ° {OxD^ + I]?+j=o ^^J-^^^y) + -^2- Comparing the 
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corresponding coefficients we have go2 = ^ ?tioo, 520 = 9ii = 5io = 5oi = 0, I2 = 1, 
P2 = 1 - .900, 7'2 = 1 - "loo 500 - goox, wliile moo satisfies m^p + mooa; = 0, and 500 
is a free parameter. Thus, we have 

A4 = (D^ + TOoo) o {D^D^ - niQoDl + goo) + (1 - 9m)Dx + Dy + 1- moogao - 9oqx , 

and the order of common obstacles is 1. 

3. Unique common obstacle and infinitely many incomplete factorizations. 
Consider factorizations of A4 of the factorization type {X){X){Y'^). We search for 
common obstacles in the form i?3 = psD^ + qsDy + r^, where p^, 93, ^3 S K . Thus, 
for some , ri3 , 03 , 63 , C3 G K we have 

Ai = (D^ + ms) o {D^ + 713) o {Dl + asD^ + bsDy + C3) + R3 ■ 

Equating the corresponding coefficients we have — —m^, 03 = 63 = C3 = 0, 
P3 = 13 — '''3 — ^, and 1713 satisfies m§ -fmsj; = 0, that is we have a unique common 
obstacle, but incomplete factorizations can be different: 

A4 = {D^ + 7713) o {D^ - 7773) o Dl + + Dy + l . 

The following lemma is used for the proof of Theorem 13.21 

Lemma 3.1 (Division lemma). LetL,AI E K[D] andSym{L) is divisible bySym{M), 
then there exist N, R E K[D] such that 

L = MoN + R , 

where either R = 0, or Sym(i?) is not divisible by Sym(M). Here R is the remain- 
der of the incomplete factorization. 

Proof. Let Sym(L) = 515*2, Sym(M) = Si. Construct a finite sequence of 77 (for 
some 77) incomplete factorizations of L of the form L = M o Ni + Qi, where 
Sym(A'i) — S2, !<«<"., and Sym((5„) is either zero or not divisible by Sym(Af). 
Start with iVi = ^2 and let Qi = L — M o Ni. If Sym(Qi) is either zero or 
not divisible by Sym(M), we stop and let N = Ni and R ~ Qi. Otherwise, 
let Ti = Sym(Qi)/Sym(M) and let N2 = Ni + T\, and let Q2 = Qi - M oT\ 
(which implies Q2 = L — M o N2). If Sym((52) is either zero or not divisible 
by Sym(Af), we stop and let N — N2, R — Q2. Otherwise, we continue in the 
same manner. Since we clearly have ord((52) < ord((5i) < ord(L), and in general, 
ord((5i_|_i) < ord(Qi), this process must stop after a finite, say 77, number of steps, 
and we have L = M o N + R, where Sym(A^) — S2 and Sym(i?) is either zero or 
is not divisible by Sym(M). □ 

Let A G K[D] and Sym(A) = 5i • ^2 • 53. It is easy to see that every com- 
mon obstacle to factorization of the type (5i)(52)(53) is the remainder for some 
incomplete factorization of the type (5i52)(53) and so it is for some incomplete 
factorization of the type (5i)(5253) (the order of the common obstacles for a fac- 
torization of the type (5i52)(53) (resp. (5i)(5253)) can be smaller that of the 
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common obstacles for a factorizatfon of the type (S'i)(S'2)(>5'3). In general, the m- 
verse statement is not true for it is more difficult to find a factorization into more 
factors. 

The following theorem states that under some conditions, common obstacles 
to factorization into two factors are the same as those into three factors. 

Theorem 3.2. Let A € K[D] and suppose Sym{A) = S1S2S3, with gcd(S'i, 6*3) = 1. 
Let U , V and W be respectively the sets of common obstacles to factorizations of 
A of the types {Si){S2S3), (S'iS'2)(<S'3), and (S'i)(S'2)(S'3). Suppose V (resp. U) is 
non-empty and the order of common obstacles in V is less than ord(S'3). Then W 
is non-empty and V = W . 

Proof. Let Ri G V he any common obstacle of type (S'iS'2)(S'3). Then we have 
ord(i?i) < ord(S'3). Let L,F e K[D] be such that 

A = LoF + Ri , (3.3) 

where Sym(L) = S1S2, Sym{F) = 5*3. Similarly, let R2 £ U, ord(i?2) < ord(53) 
be a fixed remainder with respect to an incomplete factorizations of A of type 
{Si){S2S3). Let M,G€ K[D] be such that 

A = M0G + R2 , (3.4) 

where Sym(Af) ^i, Sym(G) = S2S3. 

By Division Lemma [5t1 there exist N,Rg K[D] such that 

L = MoN + R , (3.5) 

where Sym(iV) = S'2 and Sym(i?) is either zero or is not divisible by Sym(M). 
We now claim that R = 0. Combining (|3.3I) . p.4p . (13. 5p . we have 

{MoN + R)oF + Ri = M0G + R2 , 

RoF + Ri~R2^ Mo{G-NoF) . 

Since the orders of Ri,R2 are both less than the order of F, if R were not zero, 
the symbol on the left side of the last equation would be that oi R o F, which 
would imply that Sym(i?) is divisible by Sym(M) because gcd(S'i, 53) = 1. Hence 
R ~ 0, showing that A — AI o N o F -\- Ri is an incomplete factorization of 
type {Si){S2){S3) with remainder We now show Ri is a common obstacle for 
that type. This follows easily since ii A = Mq o Nq o Fq -\- Rq is any incomplete 
factorization of that type, then A = (Mq o Nq) o Fq + i?o is one of type (5'iS'2)(-S'3) 
and hence ord(i?o) > ord(i?i). This completes the proof that V C W, which is 
thus non-empty also. If furthermore, Rq £ W, then since we have shown that 
Ri e W for any Ri E V, it follows that ord(i?o) < ord(i?i) and hence Rq is also 
of minimal order as a remainder of type (S*! S'2) (5*3), or in other words, Rq G V. 
This shows that V ^W. □ 

Example. In Theorem [321 take A to be A4 from (13. 2p . and take Si = X, S2 = X, 
and 5*3 = Y^. As we showed in the examples before Theorem 13.21 the orders of 
common obstacles of the types (5'iS'2)(S'3) and {Si){S2S3) are 1, which is less 



Refinement of Two-Factor Factorizations of LPDOs 



7 



then the order of 53. The theorem implied that the sets of common obstacles to 
factorization of the types (S'iS'2)(5'3) and (S'i)(S'2)(S'3) are the same, which accords 
with our computations in the examples before Theorem 13.21 

Example (Assumptions on the orders of common obstacles are necessary) . Consider 
operator from (g^]), where Sym(^) = X^Y^. Let Si = X^, S2 = S3 = Y . 
Then gcd(S'i,S'3) — 1. It was shown (example before Theorem 13. 2p that with 
respect to the type (S'i)(S'2S'3) = {X'^){Y'^), the operator R2 = + Dy + 1 is 
the unique common obstacle. Using similar methods, it can be shown that with 
respect to the type (S'iS'2)(S'3) = {X'^Y){Y), the operator Ri = I?-,; + 1 is a 
common obstacle. Here the hypothesis of Theorem 13.21 is not satisfied, because 
ord(ii!i) = ord(i?2) — ord(S'3) — 1. It can also be shown that with respect to the 
type (51 ) (52) (5*3) = {X'^){Y){Y), the only common obstacle is i?o = + Dy + 1. 
Clearly, i?i 7^ Rq cannot be a common obstacle of type {Si){S2){Sj,). We note 
that = Dy, i? = 1 in this example. 

Corollary 3.3. Let, in K[D\, an LPDO A have two factorizations into two factors: 

LoF = A = MoG {orFoL = A = GoM) , 

where gcd(Sym(_F), Sym(Af )) 1. Then there is a factorization of A into three 
factors: 

A^MoNoF {orA^FoNoM) 

for some N G K[D]. 

Proof. The first statement is implied from that of Theorem 13.21 For the second 
(the one which is in the brackets) we apply properties of the formal adjoints of 
LPDOs. □ 

Example (Fourth Order LPDO). Let 

L = Dl + {1 + x)DlDy + xD^Dl - x^DI ~ x^D^Dy 
+ (1 - 4a;)D^ + {x - 2x^)Dy - 2 , 

and F = Dy + x'^, M = D^ + xDy, and 

G = DlDy + D^Dl + x^D^^, + {Ax - x^)D^ + Dy - Ax^ + x^ + 2 . 

Then L o F ~ M o G, meaning that we have two different factorizations into two 
factors for the LPDO A = L o F. Moreover, one can find an LPDO N such that 
L = MoN. Explicitly, N = D^ + D^Dy - x'^D^ - 2a; + 1. Then A = M o N o F, 
meaning that A has a factorization into three factors. 

Example (Multidimensional LPDO). We have L o F = A/ o G for L = Da;Dy + 
sDj, + ts + Sx, F = Dz + b, M = D^ + t, G = D^D^ + bDy + sD^ + sb + by. It is 
also easy to see that L = M o N, where N ~ Dy + s. 
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Example (Condition gcd(Sym(i^), Sym(Af)) = 1 is necessary for Theorem 13. 3p . 
Consider L = D^Dy + j^D, + xDy + F ^ + M = + 1, 

G = D,Dy + j^^D, + ^^T^Dy - for which LoF^MoG. Here Sym(i) 

is divisible by Syni(M), but condition gcd(Sym(F), Sym(M)) = 1 fails. On the 
other hand, the Laplace invariants for LPDO L are h = —1, k = ""^fz^p^ Oj 
and, therefore, L has no factorization. 

4. Conclusions 

The main result of the paper formulated in Theorem 13.21 provides a simplification 
of the overall picture of factorization of LPDOs. 



5. Appendix 

Below is an example of how the direct approach and the approach based on The- 
orem [3]2] are different when it comes to computations. 

Let us search for factorizations of the type {X)(XY){Y) for a bivariate fourth- 
order LPDO, A — D^Dy + J2i+j=Q^ij-^x^y' ^ij ^ ^- The direct approach con- 
siders A = ]\I o N o F for some M — + m, N — D^Dy + nioDx + noiDy + noo, 
F = Dy + /, where to, nio, noi, ngo, / G K. Equating the corresponding co- 
efficients, we have 030 = 003 = 0, nio = 021 — /, noi — ai2 — to, noo — 
mf - TOa2i - /ai2 - fx - (021)3; + an, and 

= 2f.j;y - pai2 - 4,fxf + fan + '2fxa2i + fyUu - aio , 

= fy- + fa2l - 020 , 

= mail - ni^a2i - 2ma2ix - mxa2i - a2ixx + Oiia; - Ooi , 
= TOai2 -rr? - + ai2x - ao2 ■ 

= fxxy - floo + mfyai2 - rn^fy - 2fl - 2f^^f + fxxa2i + (5.2) 
+fyai2x + fxyai2 - f a2i + m'^f^ + f^m^ - f'^ai2x - 
-fa2ixx + faiix + fxttii - TO/^ai2 - 2mfa2ix + mfan - 
-fm-xa2i - 2ffxai2 - fyUix , 

An approach based on Theorem 13.31 considers A — L o F = M o G for some L = 
DlDy+Y.%,^0 hjD^Dl, F - Dy+f, M = Dx+m, G = i?.^^+E-+,=o g^jD^D^, 
where kj,f, to, gij e K . A ^ L o F imphes 030 = and /20 = 021 - /, /02 = oos, 
hi = 0-12, ho — 0,11 — ^12/ ~ 2fx, loi — ao2 ~ 003/7 ^00 = ^03/^ ~ /ao2 — ^12/2; — 
2a — 03 fy — fxx + clqi, while A — M o G implies 093 = 0, and 1720 = 0, 511 = 021, 
go2 = ai2 - m, gio = a2o, 501 = an - ma2i - a2ix, ffoo = aio - "^020 - a202;- The 
remaining conditions are 

conditions (|5.ip , 
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and two new conditions: 

= floo - maio + m'^a2o + 2ma20x - awx + mxa2o + 0202:2: , (5.3) 
= fxxy - Pam - 2/1:012/ - 'i.fxxf + /ooi - 2/^ + Jxan + (5.4) 

+ /1/O02 + fxxO-21 + ai2fxy — OOO ■ 

Thus, when algebraic manipulations only are used the difference between the 
two approaches applied to the given problem is as follows. Instead of the non-linear 
Partial Differential Equation (PDE) in two unknown variables / and m, (|5.2I) that 
we have in the first (direct) approach, the second approach implies a non-linear 
PDE in variable /, (|5.4p and another one in variable m, (|5.3p . In other words, the 
second approach gives separation of variables. 

References 

[1] I. Anderson and M. Juras. Generalized Laplace invariants and the method of Dar- 
boux. Duke J. Math., 89:351-375, 1997. 

[2] I. Anderson and N. Kamran. The variational bicomplex for hyperbolic second-order 
scalar partial differential equations in the plane. Duke J. Math., 87:265-319, 1997. 

[3] C. Athorne. A Z'^ x Toda system. Phys. Lett. A., 206:162,-166, 1995. 

[4] S.Y. Startsev A.V. Zhiber. Integrals, solutions and existence of the Laplace trans- 
formations for a linear hyperbolic system of equations. Math. Notes, 74(6):848-857, 
2003. 

[5] H. Blumberg. Uber algebratsche Eigenschaften von linearen homogenen Dtfferen- 
tialausdriicken. PhD thesis, Gottingen, 1912. 

[6] T. Cluzeau. Factorization of differential systems in characteristic p. In ISSAC '03: 
Proceedings of the 2003 international symposium on Symbolic and algebraic compu- 
tation, pages 58-65, New York, NY, USA, 2003. ACM. 

[7] T. Cluzeau and A. Quadrat. Using morphism computations for factoring and de- 
composing general linear functional systems. In Proceedings of MTNS 2006, Kyoto 
(Japan), 2006. 

[8] T. Cluzeau and A. Quadrat. Factoring and decomposing a class of linear functional 

systems. Linear Algebra and its Applications, 428(1):324 -381, 2008. 
[9] G. Darboux. Legons sur la theorie generale des surfaces et les applications 
geometriques du calcul infinitesimal, volume 2. Gauthier-Villars, 1889. 

[10] E. Shemyakova. Multiple factorizations of bivariate linear partial differential opera- 
tors. LNCS, 5743:299-309, 2009. 

[11] D. Grigoriev and F. Schwarz. Factoring and solving linear partial differential equa- 
tions. Computing, 73(2): 179-197, 2004. 

[12] D. Grigoriev and F. Schwarz. Generalized Loewy-decomposition of D-modules. In 
ISSAC '05: Proceedings of the 2005 international symposium on Symbolic and alge- 
braic computation, pages 163-170, New York, NY, USA, 2005. ACM. 

[13] D. Grigoriev and F. Schwarz. Loewy decomposition of third-order linear pde's in the 
plane. In ISSAC '08: Proceedings of the 2005 international symposium on Symbolic 
and algebraic computation, pages 277-286, New York, NY, USA, 2008. ACM. 



10 



Ekaterina Shemyakova 



[14] Z. Li, F. Schwarz, and S.P. Tsarev. Factoring zero-dimensional ideals of linear partial 
differential operators. In ISSAC '02: Proceedings of the 2002 international sympo- 
sium on Symbolic and algebraic computation, pages 168-175, New York, NY, USA, 
2002. ACM Press. 

[15] Ziming Li, Fritz Schwarz, and Serguei P. Tsarev. Factoring systems of linear PDEs 

with finite-dimensional solution spaces. J. Symb. Comput., 36(3-4) :443-471, 2003. 
[16] P.J. Cassidy and M.F. Singer. A Jordan-Holder theorem for differential algebraic 

groups, submitted to the Journal of Algebra, 1. larXivn 003.3274 April 2010. 
[17] E. Shemyakova. The parametric factorizations of second-, third- and fourth-order 

linear partial differential operators on the plane. Mathematics m Computer Science, 

l(2):225-237, 2007. 

[18] E. Shemyakova and F. Winkler. Obstacle to factorization of LPDOs. In J.-G. Dumas, 
editor, Proc. Transgressive Computing 2006, Conference in Granada Spam, pages 
435-441, Grenoble, France, 2006. Universite J. Fourier. 

[19] E. Shemyakova and F. Winkler. Obstacles to the factorization of linear partial differ- 
ential operators into several factors. Programming and Computer Software, 33(2) :67- 
73, 2007. 

[20] S.Y. Startsev. Cascade method of Laplace integration for linear hyperbolic systems 

of equations. Mathematical Notes, 83, 2008. 
[21] S. Tsarev. Factorization of linear partial differential operators and Darboux' method 

for integrating nonlinear partial differential equations. Theo. Math. Phys., 122:121- 

133, 2000. 

[22] S. P. Tsarev. An algorithm for complete enumeration of all factorizations of a linear 
ordinary differential operator. In ISSAC '96: Proceedings of the 1996 international 
symposium on Symbolic and algebraic computation, pages 226-231, New York, NY, 
USA, 1996. ACM. 

[23] S.P. Tsarev. Generalized Laplace transformations and integration of hyperbolic sys- 
tems of linear partial differential equations. In ISSAC '05: Proceedings of the 2005 in- 
ternational symposium on Symbolic and algebraic computation, pages 325-331, New 
York, NY, USA, 2005. ACM Press. 

[24] S.P. Tsarev and E. Shemyakova. Differential transformations of parabolic second- 
order operators in the plane. Proceedings Steklov Inst. Math. (Moscow), 266:219-227, 
2009. see also electronically |http: / /arxiv.org/abs/081 1. 1492 [ 

Acknowledgments. The author was supported by the Austrian Science Fund (FWF) 
under project DIFFOP, Nr. P20336-N18. 

Ekaterina Shemyakova 

Research Institute for Symbolic Computation (RISC) 

J.Kepler University 

Altenbergerstr. 69 

A-4040 Linz 

Austria 

e-mail: Ekaterina. ShemyakovaOrisc . jku. at 



